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Upon application of an external tuning param- 
eter, a magnetic state can be driven to a nor- 
mal metal state at zero temperature. This phe- 
nomenon is known as quantum criticality and 
leads to fascinating responses in thermodynamics 
and transport of the compound. In the standard 
picture, a single quantum critical point occurs at 
zero temperature, which results in a nontrivial 
critical behaviour in its vicinity. Here we show 
that in two dimensions the scenario is consider- 
ably more complex due to the enormous amount 
of quantum fluctuations. Instead of the single 
point separating the antiferromagnet from the 
normal metal, we have discovered a broad re- 
gion between these two phases where the mag- 
netic order is destroyed but certain areas of the 
Fermi surface are closed by a large gap. This gap 
reflects the formation of a novel quantum state 
characterised by a superposition of <i-wave super- 
conductivity and a quadrupole-density wave, id 
est a state in which an electron quadrupole den- 
sity spatially oscillates with a period drastically 
different from the one of the original spin-density 
wave. At moderate temperatures both orders 
co-exist at short distances but thermal fluctua- 
tions destroy the long-range order. Below a criti- 
cal temperature the fluctuations are less essential 
and superconductivity becomes stable. This new 
phenomenon may shed some light on the origin of 
the mysterious pseudogap state and of the high- 
temperature transition into the superconducting 
state in cuprates. Our results demonstrate that 
quantum phase transitions between antiferromag- 
nets and normal metals in layered materials may 
be the proper playground for search of new high 
temperature superconductors. 

Quantum criticality has become one of the most fas- 
cinating subjects in theoretical physics for the last two 
decades. At a certain point in the phase diagram, when 
the system becomes critical, quantum fluctuations are 
so strong that the metallic state is broken. Electrons 
do not behave as stable quasi-particles anymore but be- 
come heavy and short lived, hence affecting profoundly 
both thermodynamic and transport properties of those 
materials. Quantum critical points (QCP) are among 
the strongest disturbances that can be exerted upon the 
metallic state; as such they are top candidates to ex- 
plain the mysterious behavior of high-temperature su- 
perconducting cupratesi~— , heavy fermions^, or doped 
ferromagnets^. 

The beginning of the study of the quantum critical 



metals dates back to mid-seventies^. The main issue 
in dealing with quantum criticality in metals is that 
the fermions around the Fermi surface create dangerous 
massless modes which generate singularities in the com- 
putation of physical properties. The solution in these 
early works was sought in integrating the fermionic de- 
grees of freedom out of the partition function from the 
very beginning. As a result, an originating low-energy 
Landau-like quantum action describing magnetic waves 
in the normal phase of the metal — called paramagnons — 
was used to determine the upper critical dimension and 
derive renormalisation group (RG) equations. This the- 
oretical model is satisfactory for three dimensional (3D) 
systems. In two dimensions (2D), however, strong quan- 
tum fluctuations constitute a possibly major obstacle, 
capable of even invalidating the whole approach. 

An important step in the study of the quan- 
tum criticality was the formulation of a spin-fcrmion 
(SF) model designed to picture 2D metals near the 
antiferromagnetio2il£ and ferromagnetic^ instability. In 
this model, low-energy fermions interact with a collec- 
tive spin-fluctuation mode <f> that close to the magnetic 
transition is driven into criticality. 

Introducing a large number N of artificial "flavours" , 
the authors of fiefs^^ developed a technique similar 
to the one well known in the context of strongly cou- 
pled superconductivity. It consists of renormalising all 
the fermion propagators while neglecting the vertex cor- 
rections. Proceeding in this way, one reproduces in 
the limit N — > oo the results of the Hertz theory^, 
though the physical picture is enriched by intriguing 
renormalisations^— of Fermi velocities and vertices in 
2D for finite N. As such, for a quite long time of more 
than ten years, the (anti)ferromagnet-normal metal tran- 
sition was believed to have been understood with the 
critical behavior characterised by scaling relations that 
resemble those for second-order phase transitions. 

However, as recently unveiled^ 2 - - — , this ansatz is in- 
complete as the conventional 1/Af-expansion unexpect- 
edly misses contributions in both perturbation theory 
and renormalisation group that actually are not sup- 
pressed by 1/N — thus invalidating the picture. As a 
result, with no controlled expansion at hand, researchers 
could not come to a definite conclusion about the char- 
acter of the quantum phase transition and the problem 
has remained open. 
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Model and method of calculations 

In this Article, we revisit the issue of quantum 
antiferromagnet-normal metal transitions in 2D models 
of itinerant electrons from the perspective of a novel gen- 
eral theory that shall lead us to the conclusion that the 
physics of the transition is considerably richer and more 
interesting than it has been thought so far. We demon- 
strate that within a slightly modified version of the SF 
model of Refsj*^, the coupling of the bosonic spin mode 
to the electronic spins generates at the QCP a pseudogap 
in the spectrum that corresponds to an order completely 
different from the original spin-density wave (SDW) . This 
new state may be understood as a superposition of d- 
wave superconductivity and an electronic quadrupolc- 
density wave and its emergence around the QCP consti- 
tutes an unexpected outcome of our theory. Interestingly, 
our formalism shows a certain analogy with the theory of 
Anderson localisation by disorder—. The structure of the 
both theories relies on a summation of ladder diagrams 
of crossing and non-crossing subtypes, that characterise 
the emerging effective collective modes and their inter- 
action. In both cases, the low energy physics is finally 
captured in terms of a non-linear a-model. Below, we 
present a sketch of the derivation and refer the reader to 
the supplementary material for details. 

Within the SF model, the physics of electrons inter- 
acting via critical bosonic modes is described*^ by the 
Lagrangian L = + with 

L. 4 , = t/j* [d r + e (-»V) + AH V, (1) 

Lt=\<t>D-*<t>+ 9 - (0 2 ) 2 . (2) 

Herein, is the Lagrangian of fermions that propagate 
in the fluctuating field <p representing the bosonic spin 
excitations modeled by the Lagrangian L^. The elec- 
tronic spectrum e(p) in is assumed to lead to a Fermi 
surface like the one depicted in Fig. [T] The Lagrangian 
L<f, is a quantum version of the Landau expansion in the 
vicinity of a phase transition. 

We define the spin-wave boson mode D~ x entering 
Eq. §2$ around the QCP through its Fourier transform 

D- 1 (u,q)=LU 2 /v 2 s + (q-Qj 2 + a (3) 

where v s is the spin-wave velocity, a is a "mass" charac- 
terising the distance to the QCP [At QCP a — 0, while 
a > on the metallic side.], and Q is the ordering wave 
vector in the SDW phase. Keeping in mind possible ap- 
plications of the model to high-T c cupratcs, we choose 
Q = (tt.tx). The points on the Fermi surface connected 
by the vector Q are the hot spots in the model and, close 
to criticality, the most interesting physics is formed in 
their vicinity. Figure Q] illustrates that there are eight 
hot spots on the Fermi surface. It is implied that the 
coupling constant A is small in the sense of A 2 <C vp , 
where v is the Fermi velocity and po the radius of curva- 
ture at the hot spots (see Supplementary Information for 




FIG. 1: Brillouin zone and Fermi surface for the spin- 
fermion model. Electrons at hot spots connected by the 
vector Q interact via a critical bosonic mode. The vec- 
tors Qi,2 modulate the amplitudes of particle-hole pairings. 
Inset: Definition of the angle S that controls the theoretical 
approach. 

the details). Also, the quartic <^ 4 -term in is usually 
neglected for a > 0. 

Following Refs^ i 10 ' 14 , we might assume that the 
fermionic field t/j had 2N components, where 2 arises due 
to spin and N is the number of artificial fcrmion flavours. 
However, in accord with the conclusion of Reb 1 ^, the 
large- N limit does not help to justify approximations and 
control the theory and, therefore, we shall take here the 
physically relevant value N = 1. Instead, we keep the 
theory under control assuming that the Fermi surface 
has such a shape that the Fermi velocities v 12 of two 
hot spots connected by vector Q are close to being par- 
allel to each other, 

S < 1, (4) 

with the angle 5 defined in the inset of Fig. [1] Even 
though the constraint (j4|) itself is in principle not unre- 
alistic, indicating the tendency to a nesting of the type 
represented in Fig. [TJ we believe that the results of the 
present analysis can be applicable for an arbitrary shape 
of the Fermi surface at least qualitatively. 

The Landau damping modifies^ the form of -D -1 (w, q), 
Eq. ([3]), adding to the latter the term j\ui\ with 7 = 
(2A) 2 /(7rv 2 sin5) and v = |vi j2 |. In the limit g}, the 
Landau damping is strong, leading to a "weak coupling" 
limit of our theory, and dominates over the w 2 -term in 
the bare propagator D^ 1 . 

In the spirit of the approach of Ref we first integrate 
the partition function Z = J exp{— J L}D(f>Dip over the 
field <f> neglecting the quartic term in L$. As a result, we 
obtain a model of electrons with an interaction described 
by the function D (w,q) including the Landau damping 
term. Next, we single out those slow pairs that corre- 
spond to the mean field order parameters and derive the 
mean field equations. The order parameter found from 
these equations is strongly degenerate against SU(2) ro- 
tations, which gives rise to gaplcss excitations that finally 
are effectively described in terms of a non-linear er-modcl. 
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FIG. 2: Chequerboard structure. The quadrupole density 
amplitude (normalised to values between ±1) is represented 
in real space. It is incommensurate with the square atomic 
lattice of the compound. 



The properties of the latter are studied using the renor- 
malisation group (RG) technique. 



Mean field equations and pseudogap state 

The mean field approximation leads to a superposition 
of particle-particle and particle-hole pairings, 



,,]>]> 



( ia 2)al3 ^a,p1pp,-p) + tip (<*a/jV'a J P^S,-p) i ( 5 ) 



with the momentum p located at hot spots opposite 
to each other on the Fermi surface. In Eq. ([5]), 02 is 
the Pauli matrix for the electron spin and C pp (eg) the 
amplitude of the particle-particle (particle-hole) pairing. 
The pairings of the type in Eq.([5]) are purely singlet and 
thus do not lead to any spin order. The signs of coef- 
ficients Cp P//ph at neighbouring hot spots on each con- 
nected piece of the Fermi surface are opposite, indicating 
a d-wave-like structure of the gap in the electron spec- 
trum. Therefore, neither local charge nor current den- 
sity modulations arise. However, as a consequence of 
the electron-hole pairing, the rotational symmetry of the 
electron gas is broken and a state with a finite coordinate- 
dependent quadrupole density is formed. The electronic 
quadrupole density, which is proportional to the pair- 
ing amplitude, oscillates with the vectors Qi and Q2, 
sec Fig. [TJ connecting two hot spots at ±p. These vec- 
tors Qi_2 are considerably smaller than the SDW order- 
ing wave vector Q and a chequerboard structure as de- 
picted in Fig. [3] forms (c.f. RefJ^). We dub this new 
structure quadrupole- density wave (QDW). At the QCP 
the parameter a vanishes, a = 0, and the propagator 
D(oj,q), Eq. ((3]), diverges at q > Q and ui — > 0. The 
latter ultimately allows in 2D for the electron-hole pairing 
in the case of a generic Fermi surface without a nesting. 
This pairing survives also for finite but small a. 

The general solution O of the mean field equations for 
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FIG. 3: Gap function. The mean field solution b(e, T) as 
a function of the Matsubara frequency e and temperature T. 
All energies are measured in units of V. 



the order parameter at a given hot spot may be repre- 
sented in the form 0(e) = b(e)u with u being an arbitrary 
SU(2) unitary matrix, u + u = 1, dct u = 1, and 6(e) a real 
positive function of the fermionic Matsubara frequency e. 
After a rescaling e — > eT, b — > bT, and T — > TT, where 
r = (3A/8) 7rsin<5, one obtains at criticality (a = 0) 
a set of remarkably universal self-consistency equations 
independent of the parameters of the model, 



b(e)= 



He) 



cos9(e') 
y/to (e -e<y 
sin9(e') 



t\,-\ o • 2 /9(e + o>) - 9(e) 

SI (u>) = 2nT sin ' 



(0) 



-1/2 



where sin 9(e) = /(e) [b 2 (e) + / 2 (e)] . The functions 
6(e) and /(e) are_by construction even, 6(e) = 6(— e), and 
odd, /(e) = — /(— e), respectively and Q is a rescaled 
bosonic Matsubara frequency. Note that the function 
/(e) replaces the frequency term e in the bare fermion 
propagator. 

A quick glance at Eqs. © reveals the trivial solu- 
tion 6(e) = 0, leading to Q(p) = \u>\ and /(e) = 
sign(e)(|e| + -■^/f¥f). This solution is well known as it 
corresponds to the one-loop self-energy correction o 9 ' 10 ! 17 
to the bosonic and fermionic propagators. Here, of a 
greater interest is the existence of a nontrivial so far 
unanticipated energy-dependent solution 6(e). It can be 
computed numerically and its dependence on e and T is 
shown in Fig. [3] We have checked that the free energy 
corresponding to the nontrivial solution is lower than the 
one in the case of the trivial scenario with 6(e) = 0. One 
should keep in mind, though, that the dependence of the 
gap on real frequencies does not immediately follow from 
the solution 6(e) but should be found from an analytical 
continuation. 

The characteristic value of 6(e) is of order T, im- 
plying that it scales linearly with the interaction con- 



FIG. 4: Gap as a function of the position on the Fermi 
surface. The gap 6(0, p) is essentially non-zero only in the 
vicinity of hot spots. The order parameter has opposite signs 
at the hot spots located on the same pieces of the Fermi sur- 
face, which corresponds to a d-wave-like symmetry. 

stant A 2 . This is in a sharp contrast with the exponen- 
tially small values of the gap encountered in conventional 
superconductors^. Therefore, one can expect at QCP 
much higher values of the gap than those obtained for 
non-singular interaction functions used in BCS theory. 
Eqs. ([6]) have been obtained linearising the electron spec- 
trum near the Fermi surface and their solution formally 
does not depend on the position on the Fermi surface. 
In fact, the order parameter O depends not only on the 
frequency e but also on the distance from the hot spots, 
decaying at momenta of order j/v. The length of the arc 
of the Fermi surface under the gap, however, should be 
calculated taking into account the curvature of the Fermi 
surface at the hot spots. As a result, the gap is finite only 
in the vicinity of the hot spots, while the "cold regions" 
of the Fermi surface remain gapless and their contribu- 
tion to physical quantities is typical for the conventional 
Fermi liquid. Schematically, 6(0, p) on the Fermi surface 
is depicted in Fig. [4] 

The resulting momentum dependency of &(0, p) on the 
position of the Fermi surface can in principle be de- 
tectable in STM or ARPES measurements. Actually, as 
the function b depends on the Matsubara frequency and 
the position on the Fermi surface, it describes a pseudo- 
gap rather than a real gap in the spectrum. The corre- 
sponding state obtained from the non-trivial solution of 
Eqs. ([6|) can be called pseudogap state. 

The matrix u reflects degeneracy of the order parame- 
ter and can be parametrised as 

U =( A A* + A t ) with |A+| 2 + |A-| 2 = 1. (7) 

The complex numbers A + and A_ should be interpreted 
as order parameters for the superconducting and QDW 
order, respectively. In contrast to the conventional su- 
perconductivity where electron-electron pairs are formed, 
we have here quartets consisting of two particles and two 
holes, see Fig. [5] Depending on the relation between the 
horizontal and vertical coupling one of the pairings is 
more favourable but, if they are equal or one considers 
fluctuations, the entire quartet should be dealt with. 

The nature of the particle-hole pairing in our theory is 
different from those conjectured in SU(2) theories on the 
basis of symmetries of t — J models^. The QDW order 
found here arises without generating any modulations of 



FIG. 5: Quartets formed in the pseudogap state. Pair- 
ing types of electrons and holes at opposite hot spots for (a) 
quadrupole-density wave (QDW) order and (b) d-wave super- 
conductivity. 

the local charge density or currents and thus cannot be 
directly identified in STM or ARPES measurements. At 
the same time, charge density modulations can be in- 
duced by adding inhomogeneities into the system. 

Fluctuations and phase transition into 
superconducting state 

The degeneracy with respect to rotations of the ma- 
trix u leads to gapless excitations. They destroy the 
long-range order and smear the transitions. This phe- 
nomenon is well-known in the context of the condensa- 
tion of a Higgs boson and of the description of the fluc- 
tuations around ground states. Following a similar route 
here but taking into account the more complex symme- 
tries of the model, we study the contributions of these 
excitations to thermodynamics by deriving the proper 
non-linear cr-modcl in a way developed in the localisa- 
tion theorjiLii. At T = 0, the model is 2 + 1 dimensional 
and the contributions of fluctuations of u as well as of 
the fluctuations of the amplitude of the order parameter 
are converging and small in the limit ((3]). As a result, 
the mean field equations ([6]) are applicable. Fluctuations 
at finite temperatures are more dangerous because the 
effective dimension is reduced to d — 2. The derivation 
of the effective free energy functional F [u] that describes 
the fluctuations yields 

Flu] It' 

= - I tr [Vu+Vm + k 2 u+t 3 ut 3 ] dR . (8) 

In Eq. flgj), t = c(T)T/\ 2 , where c(T) is a monotonous 
function of the order 1 at low temperatures. The latter 
vanishes at the mean field transition temperature To ~ T 
of the pseudogap state. Further, k ~ ~f 2 v/m with m be- 
ing the electron mass, and R = {xj sin(<5/2), yj cos((5/2)} 
where x and y are the original coordinates directed along 
the diagonals of the Brillouin zone in Fig. [T] Note that 
the coupling constant t does not contain sin S and is small 
for all temperatures that are below but not too close 
to Tq. The matrix T3 stands for the Pauli matrix in the 
space of the matrix u, Eq. ([7]). Its presence in the a- 
model breaks the symmetry between the superconducting 
and QDW states. As a result, a superconducting order 
with a large gap ~ Y is more favourable at the minimum 
of F[u]. 
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It is well known that fluctuations in two-dimensional 
non-linear er-modcl at finite t and k — produce log- 
arithmically divergent contributions destroying the long 
range order. This means that in our case they destroy 
both the superconducting and QDW orders. A finite k 
serves as the infrared cutoff in the logarithms, so that 
the superconducting order can be stabilised below a crit- 
ical temperature T c . The RG method is a standard tool 
for studying properties of the cr-modeli^. Integrating out 
step-by-step all momenta exceeding k in the first loop 
approximation we come to the following expressions for 
the effective coupling constants £(ko) and k(kq) 

t( K ) = t (l-^L\n-) \ (9) 
V l07T vk J 

- 2 M^(l-^hr-) , (10) 

where to and kq are the bare values written below 
Eq. (fl~5| and T plays the role of the upper energy cutoff 
in the cr-model. 

Formally, Eq. © is applicable as long as t(n) <C 1 
where superconductivity is stabilised. Upon increasing 
the temperature, the bare coupling t and, hence, t(n) 
grow, while Kq decays. These tendencies imply that fluc- 
tuations of the order parameter u become strong but, at 
the same time, the anisotropy between the superconduct- 
ing and QDW vanishes. This should lead for t(n) ~ 1 to 
a phase transition into a "disordered" (pseudogap) phase 
where no long-range order exists anymore but the pseu- 
dogap b(e) is still finite and remains large up to consid- 
erably high temperatures. The pseudogap state should 
exhibit combined properties of d-wave superconductiv- 
ity and QDW although sharp dependencies are neces- 
sarily smeared by thermal fluctuations. An estimate for 
the transition temperature T c between the superconduct- 
ing and pseudogap phase is provided by the equation 
= 1. It is a question of future studies how various 
features of the pseudogap state can be probed experi- 
mentally, whereas the pseudogap itself can certainly be 
seen using STM and ARPES techniques. It is also evi- 
dent that there remain strong superconducting fluctua- 
tions above T c . Note that in contrast to the BCS theory 
of superconductivity, T c is related to the value of the gap 
r in quite a nontrivial way. 

Phase diagram 

We summarise our findings with the help of the phase 
diagram in Fig. [6j The red dashed line denotes the 
antifcrromagnet-normal metal phase transition in the ab- 
sence of interaction between the spin-wave modes and 
electron spins. It is the region to its right that we have 
studied in the present Article. We identify the pseudogap 
region and represent the crossover to the normal metal 
state by the blue line. The green superconducting re- 
gion is a part of a more broad pseudogap region. Moving 




bare QCP 



FIG. 6: Phase diagram for the spin-fermion model. In 

this picture, AF denotes the antiferromagnetic (SDW) state, 
SC is the phase of the d-wave superconductivity, and PG the 
pseudogap state. The dashed line represents the solution of 
the equation a(T) — 0. The question mark "?" indicates that 
the present consideration is not sufficient to identify the phase 
in the region between AF and SC. 

to the right (increasing the bosonic mass a in Eq. ©) 
makes the spin-fermion interaction less singular and ef- 
fects of curvature of the Fermi surface more pronounced. 
As a result, QDW pairing is suppressed but one still ob- 
tains the d-wave superconductivity, although with a lower 
transition temperature. At large a, the pseudogap and 
superconducting states merge. 

Although our results cannot be used to the left of the 
dashed line (a < 0), by assumption, the antiferromag- 
netic region exists within the SF model, Eqs. (|1I2[) , for 
not very small \a\. At the same time, the pseudogap re- 
mains finite for a while after crossing the dashed line, 
because it would cost a large energy for it to vanish. 
Therefore, there must be a region of a < in which the 
pseudogap is still finite and this region extends to the 
left until the phase transition (crossover) to the antifer- 
romagnetic state sets in. At low temperatures, there is 
the possibility for an interesting scenario that, when go- 
ing to the left, the coefficient k 2 , Eq. ((8|), changes its 
sign. This would correspond to a transition from the 
superconductivity to QDW followed by a transition into 
the antiferromagnetic state. At the same time, we cannot 
exclude a direct transition from the superconductivity to 
the antifcrromagnet and the situation will only be clari- 
fied by an explicit study of this region. We have marked 
this region by the question mark in Fig. [6l 

Remarkably, the phase diagram represented in Fig. [6] 
looks similar to that of the cuprates (see, e.g. Refs . 20 ' 21 ). 
The large values of the pseudogap can clearly be under- 
stood from the linear dependence of the pseudogap on 
the effective electron-electron interaction obtained in the 
present work. Of course, in order to make a definite con- 
clusion one should perform a detailed comparison of the 
theoretical prediction made on the basis of the SF model 
and experimental data. However, our results are already 
consistent with important experimental evidence. For in- 
stance, the d-wave like shape of the pseudogap has been 
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revealed by ARPES measurements 2 ^, while STM images 
reveal the breaking of the rotational symmetry of the 
electronic structure 2 ^. 

There are other interesting phenomena observed in 
cuprates, such as linear-T resistivit y 22 ' 23 or electron den- 
sity modulatio n 24 ' 25 . They do not immediately follow 
from the SF model, Eq. (J]) and ©, and apparently the 
model should be extended in order to gain the capacity to 
explain these effects. We think that the presence of inho- 
mogeneities well known in cuprates^ 5 - can be important 
for building a more complete theory. At the same time, 
we believe that these additional effects are not directly 
related to the origin of high-T c superconductivity. 

Outlook 

We have shown through a controlled solution of 
the spin-fcrmion model that the antiferromagnet-normal 
metal quantum phase transition in two dimensions is a 
complex phenomenon belonging to an unique universal- 
ity class different from those well known in the theory 
of phase transitions. Our most remarkable finding is the 
emergence of a pseudogap state at and near the critical 
line. It is characterised by a 2 x 2 matrix order parameter 
that physically corresponds to a superposition of super- 



conducting and quadrupole density wave orders. There- 
fore, modulations of the local electron density and cur- 
rents are absent, whereas the rotational and translational 
symmetries are broken. The parts of the Fermi surface in 
the vicinity of the hot spots are closed by the large gap 
and its dependence on the position on the Fermi surface 
resembles the one of a d-wave. At the same time, a long- 
range order is destroyed at moderate temperatures. At 
low temperatures, the system becomes a d-wave super- 
conductor with a large gap. 

A rough comparison of our findings with several robust 
features of the high-T c cuprates suggests that the theory 
presented here may be a good candidate for the explana- 
tion of phenomena related to high temperature supercon- 
ductivity. We emphasize that the model presented here 
can be fully studied in a controlled way within the frame- 
work of quantum field theory, a feature quite unique in 
the field of cuprate superconductivity and quantum crit- 
icality. Moreover, the spin-fermion model can also ef- 
ficiently be studied numerically using the Monte-Carlo 
method due to the possibility to avoid the sign problem 
here 2 ^. All this allows us to expect in the nearest fu- 
ture a considerable progress in understanding physics of 
cuprates and other materials exhibiting quantum phase 
transitions. 
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I. SYMMETRIES OF THE SPIN-FERMION MODEL. 



In the earlier considerations 8 ' 10,14 , the spin-fermion model was analysed assuming an arbitrary number N of artificial 
electron species or "flavours" . The general idea was that the limit of large N should help to justify diagrammatic 
expansions. Only recently, contributions that are, contrarily to what one would expect from conventional 1/N- 
cxpansions, not small in the limit TV — » oo have been found 14 . In this work, we confirm the conclusion that the 
large- -/V limit does not help to solve the model. Instead, it turns out that one can use in the calculations a different 
small parameter to control the theory and can thus safely put N = 1 corresponding to physical reality. Nevertheless, 
we keep an arbitrary TV in this Supplementary Information (SI) in order to demonstrate explicitly what happens as 
TV becomes large. 

The Lagrangian L of the spin-fermion model has been written in Rcf. 1 in a form that allows one to see explicitly 
symmetries of the system and we adopt it here, 



L = Lf + Lif> , 



i/i=i 

N 
~2 



L 



]T xL + (fir ~ ^V) xl + A (x[ + fcx l 2 + X l 2 + fcx[) 

a=l 

(?) 



(11) 
(12) 

(13) 



The Lagrangian (|11H13|) describes the fermions near the hot spots of the Fermi surface, see Fig. 1 of the main text, 
that interact via a collective bosonic mode. The 2iV-component anticommuting vectors Xa represent the electron 
fields with 2 being due to spin and N the number of the flavours, and the vectors v l a are Fermi velocities at the hot 
spots. Herein, the superscript I numerates the pairs of the hot spots connected by the vector Q ~(tt,tt) of the spin- 
density wave, while the subscript a distinguishes the single hot spots inside each pair. We choose coordinates along 
the diagonals of the Brillouin zone. The three-component vector field <f> describes the bosonic spin-wave modes with 
velocity v s that couple to the electronic spin, represented in Eq. (|12p by the operator a which is the vector of Pauli 
matrices. The parameter a(T), the bosonic "mass" , determines the distance from the critical line in the phase diagram 
and vanishes at criticality. In the spin-fermion model, we consider fields x an d <P depending on (1 + 2)-dimensional 
coordinates X = (r, r) combining imaginary time r, which varies between and the inverse temperature ft = 1/T, 
and coordinates r for the effectively two spatial degrees of freedom. The parameter a(T) determines the distance from 
the critical line of the phase diagram. 

The spectrum of the fermions written in Lagrangian Lf, Eq. (|lip . has been linearised. From the Fermi velocities 



= (V X Vy) , V2 1 = (-V X Vy) , 



(14) 



in the hot-spot pair I = 1, we obtain all other by inversion and rotation by 7r/2. The linearised form is not always 
sufficient for the computation of physical quantities and a slightly more general Lagrangian taking into account 
curvature effects may become necessary. We shall include the curvature corrections shortly. 
For the purpose of clear and compact notations, we introduce 32/V-component fields >£, 



Herein, the label r stands for the particle-hole space and x is a 167V-component field, 



, x + = (x 1+ x 2+ ) } 



(15) 



(16) 



The superscript L = 1, 2 distinguishes the pairs of hot spots on the Fermi surface that are connected by Q (for L = 1) 
from those pairs connected by Q' = (— ir,ir) (for L = 2), cf. Fig. 1. Finally, the components of 



x\ 

Xi L 

V V xi L / E 



x L+ = 



ix[ + xi 



h W + xf) s )^ 



(17) 
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are the spinor fields Xa appearing in the original formulation (|12p . In these equations, * maps the anticommuting 
variable x onto its conjugate and t denotes transposition. We use the convention that (?7i%)* = ViVj an d (vt)* = ~Vi 
for arbitrary anticommuting variables rji. The symbols £ and A refer to the pseudo-spin spaces corresponding to the 
subscript a (distinguishing hot spots within one critical pair) and the superscript / (distinguishing the two pairs of 
hot spots) , respectively. Thus, for N = 1, the field ^ is a tensor in the 32-dimensional space T<8>L<g)A(g)E(g>(7 with a 
being the physical spin. In the following, we shall use besides er the Pauli matrices f, A, and X in the corresponding 
factor spaces. At the same time, all relevant physical quantities are diagonal in L-space. 

Besides the hermitian conjugate \E' + of the vector it convenient to introduce in addition a "charge" conjugate 
field ^ defined as 

<F = (Ctf)* with C = y ^ 1<J Q 2 ^ = -r 2 a 2 . (18) 

The matrix C satisfies the relations C*C = 1 and C = C*. It is clear that 

* = f + T 3 . (19) 
The notion of charge conjugation is naturally extended to arbitrary matrices M (X, X') as 

(V(X)M(X,X')-$(X')) = -($(X')M(X,X')$(X)) . (20) 

It is easy to see that 

M (X, X') = CM 1 (X', X) C* = CM T {X, X') C* . (21) 
We shall refer to matrices satisfying the relation 

M = —M (22) 

as anti-selfconjugated. Such matrices frequently appear when dealing with quadratic forms for anti-commuting vectors. 

All the notations (|15H21[) are very similar to those used in Ref 16 . Using the multicomponent field VP and including 
also the curvature corrections into the spectrum of the fermions, we write the action Sf that corresponds to the 
Lagrangian Lf, Eq. (fT2|) . as 

S f m = J * (X) (n + n cn „ + AEict*^ (X)j * (X) dX (23) 



with Ho = —d T + iVV and 'Hcurv = — 1~ 3 



2 

V x V 



(24) 



Herein, m is the effective fermion mass at the hot spots and the velocity operator V, written as a matrix acting in 
the space of x, Eq. (fTo) . has the form 

(v + +£ 3 v_)A 3 \ 

V -l (v + +S 3 v_)A 3 L {2b > 



with 



T l=l | „i = l v /=l _ yl=l 



v+ = 1 2 Z =(0,v y ), v_= 1 2 4 =(v x ,0) (26) 

1=2 , 1=2 1=2 _ 1=2 

v+ = 1 2 2 = (-^,0), v_ = 1 2 2 = (0,v x ) . (27) 

In Jicurv only the momentum component directed along the Fermi surface is kept, which is sufficient in the weak 
coupling limit 

A 2 < vp , (28) 

where v and po are respectively, the modulus of the velocity and the radius of the curvature of the Fermi surface at 
the hot spots. The velocities v± are mapped onto v± by rotation over the angle ir/2. 
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All parts of the Hamiltonian, Ho, H CU i V , and A£i<7*</> (X), are anti-selfconjugate, an immediate consequence of 

B T = -d T , V = -V , and a 1 = -a 1 for i = 1, 2, 3 . 

Generally, the action is invariant under transformations 

* -> [/* (29) 

where f is a 32Af x 32Af unitary matrix with unity blocks in spin and particle-hole spaces, while the diagonal blocks 
in the L (g) A <S> E-space, the space in which the matrix V, Eq. (|25|) . acts nontrivially, are arbitrary. In other words, U 
consists of eight independent U(AT) matrices, reflecting the independent particle number conservation at each single 
hot spot. If we neglect the term H curv containing t 3 in the action Sf, the latter becomes invariant under arbitrary 
particle- hole rotations, thus enlarging the symmetry group for transformations of the type (|2"5|) to U(2iV). 

This \J(2N) symmetry of the linearised action is crucial for physical properties of the model and we shall see that 
its spontaneous breaking in the vicinity of the QCP leads to a rather intricate order parameter. 



II. MEAN FIELD APPROXIMATION 



A. Self-consistency equation 

We computate the partition function 



exp 



-S f [*,</>]- / L^} dX 



D^D(j) (30) 



first by integrating over the bosonic field 4>. This is possible in a closed form only if we neglect the quartic interaction 
in L(j), Eq. (|13|). which constitutes our basic approximation. In conventional theories of phase transitions it is not 
legitimate to neglect the quartic term since perturbation theory in g breaks down at the critical point at and below the 
upper critical dimension. In the present model, one encounters logarithmic divergencies in the perturbation theory 
at zero temperature 8 . However, renormalisation group (RG) studies of the <j> theory show a logarithmic decay of the 
coupling constant g and therefore following Ref. 8 the quartic term is usually neglected. 

The situation is more tricky at finite temperatures because the </> 4 model does not have a phase transition in 
two dimensions, while neglecting the quartic term the system would undergo the phase transition. This problem is 
discussed later with a conclusion that the quartic term can still be neglected at the cost of assuming that effectively 
the parameter a (T) should always remain finite at finite T but vanishing in the limit T — > 0. Such a behaviour of 
a (T) should mimic the smearing of the phase transition in the absence of the quartic term. 

Then, integrating out the bosonic field <fi yields 

Z = j 'cxp (-5 [*])£>*, S = So + S curv + S int , (31) 

where the noninteracting terms for the linearised kinetic spectrum and curvature corrections read 

So (X) Ho* (X) dX , S CUIV [*]=[* (X) H CU „V (X) dX . (32) 

The fermion-fcrmion interaction 

A 2 f 

Smt [*] = -y / (* (X) Ekt** (X)) D(X- X') (* (X r ) Eia** (X')) dXdX' (33) 
includes as interaction potential the original bare bosonic propagator 

D(X-X') = T W exp (-iw (r - r') + zq(r - r')) D («, q) ( 34 ) 

with D _1 (cj,q) = N(uj 2 /v 2 s +q 2 + a) . (35) 

lu = 27rTm, m = 0, ±1, ±2, . . ., are Matsubara bosonic frequencies. Note that in SI we measure the electron momenta 
from the hot spots. As a result, the propagator D (w,q) in Eq. (I3"5"j) has formally been shifted by the vector Q with 
respect to the one in Eq. (3) of the Article. 
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(a) 



-» — 



(b) 



» i » > i » 



FIG. 7: Particle-particle and particle-hole non-crossing ladder diagrams. 
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FIG. 8: Renormalisation of the spin-wave propagator by particle-hole bubbles. 



In principle, one can study the model denned by Eqs. (|31H35j) using an expansion in S[ nt . The bare Green function Go 
for the Hamiltonian "H used in this type of perturbation theory is written in Fourier space as 



G 1 ( £ >P) = i£ 



(36) 



with e = 7r (2n + 1) T, n = 0, ±1, ±2, . . ., denoting a fermionic Matsubara frequency. 

This approach has been used in the previous publications 9,10,14 . It has been found in Ref. 14 that anomalous (in 
terms of the expansion in 1/N) contributions come from particle-particle and particle- hole ladder diagrams of the 
type represented in Fig. [7] 

These diagrams describe superconducting and some kind of insulating fluctuations demonstrating a tendency to a 
corresponding particle-particle and particle-hole pairing. Formally, these diagrams resemble "cooperons" and "diffu- 
sons" appearing in the localisation theory 16 . More complicated diagrams correspond to an interaction between these 
effective modes. 

A very efficient way to sum the contributions of all these diagrams is to derive a non-linear a-model 16 and study 
fluctuations of an "order parameter" using this effective field theory. 

Below, we follow a similar route. We investigate the model by first writing self-consistent mean field equations, then 
solving them, and finally studying fluctuations. Of course, one can describe the fluctuations within a perturbation 
scheme again. However, now the expansion will be performed near another minimum. This is a standard situation in 
models where a symmetry of the original Hamiltonian is broken in a certain region of parameters. 

The first step consists of replacing the 'J -interaction by a quadratic term with coefficients to be determined 
in a self-consistent way: The effective action for the fermion-fermion interaction shall contain a bosonic propaga- 
tor D c g (X — X') — physically the spin susceptibility — that is renormalised by particle-hole bubbles, see Fig. 
that in turn are composed by renormalised fermions. In order words, we develop a self-consistency scheme for both 
fermions and bosons. 

The renormalised propagator D e g (uj, q) can be written as 



D-£ (u, q) = D- 1 (u, q) - n ( W> q) , 



(37) 



n( w ,q) = -A 2 T^ yTr[G(e,p)S 1( T t G(e + W ,p + q)S 1( T t ] A (38) 
In Eq. (|38p. G (e, p) is the Green function of the renormalised fermions and Tr is the trace of the 32iV x 32iV matrices. 
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In order to write equations for G (e, p), we replace the interaction term [$] by a quadratic form, 

S int [¥] = y /" £> eff (X - X') Tr [* (X') * (X) Sit?'* (X) * (X') £ 1( r 4 ] dXdX' (39) 

~ i y Tr [* (X') * (X) M {X, X')] dXdX' 

+ (2A)~ 2 y Z^ 1 (X - X') Tr [M (X, X') Eict*M (X 'AT) Sict 4 ] dXdAT' . (40) 

In the spirit of the mean field approximation, the matrix M(X, X') should be found in a self-consistent way from the 
equation 

M (X, X') = -2i\ 2 D cff (X - X 1 ) (E 1( i 4 * (X) $ (X 1 ) £i<r 4 ) off (41) 
where the symbol (. . .) s is the mean field average, 

oxp(-S off M)£>tt 

\- /cff Jexp(-S eS [^})D^ ■ 1 J 

Herein, the effective action S e g [VP] reads 

5* cff [*] = y * (X) Heff * (X) dX, (43) 

H cff = -9 T + iVV - iM (X) 
with the operator M (X) acting on the field \I> (X) as 

M (X) = [ M (X, X') * (X') dX' . 



The factor of 2 in Eq. (|4"Tj) arises due to the two possible pairings of "Fock-" and "Cooper" -type that are in fact equal 
to each other. 

Now, we are in a position to define the Green function G(X,X') as the propagator for the action S e s, Eq. (|42j) . 

G(X,X') = -2(*(X)*(X')) gS . (44) 

It satisfies the equation 

(d T - iVW + iM^j G (X, X') = 5(X-X') . (45) 

Using G(X,X'), the self-consistency equation (|4"Tj) takes the form 

M (X, X 1 ) = i\ 2 D cS (X - X') Ei,r 4 G (X, X 1 ) H^ 1 . (46) 



Equation (|46| is the basic equation that determines all the unusual physics at and near QCP. It is an analog of 
the BCS equation in the conventional theory of superconductivity 18 . However, the singular form of the interaction 
D e ff (X — X') leads to a more complex mean field solution and, as a result, to richer physics than that given by the 
BCS theory. 

B. Solution of the self-consistency equation 

By Eq. ([41j) . the matrix M (X, X') is anti-sclfconjugate (with the variables X X' included in the transposition), 

M (X, X') = —M (X, X 1 ) . (47) 
We look for the solution M (X, X') to Eqs. (|45I46|) using the ansatz 

M(X,X')=M (X-X') , (48) 
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+ /V"\ + 

'c' c -c* C ' c' c' c ■ c c' 



FIG. 9: Summation of the non-intersecting diagrams for the self-energy of the fermionic Green functions leads to the self- 
consistency equation. 



which is Fourier transformed as 

Mi 



(X - X') = [ M (e, p) e Hr-T')-i P (r-r')JP . (49) 

The self-consistency equation (|46|) thus becomes 



M (p,s) = *A 2 TV f D off (e-e / ,p-p / )E 1 ^G(e / > p / )Si5*^j (50) 
J (2tt) 



G(e,p) = (ie- Vp + iM (e,p)) ' . (51) 



where 



This Green function should also be used in Eqs. (|38|) for the propagator D c g so that the system of the self-consistency 
equations is closed. Note that 

M (e,p) = -M (-£,-p) . (52) 



as a consequence of the anti-selfconjugation property of M, Eq. (|47|). 

Eqs. ([50l l5Tj) can also be obtained by summation of all non-intersecting diagrams for the self energy like those 
depicted in Fig. 

The presence of the spin matrices a in Eq. (|40[) prevents the solution from being invariant with respect to rotations 
in the spin space. Parametrising the spin structure of Mo(e, p) as 

M (e, p) = Q (e, p) + ia'S (e, p) , (53) 

where Q(e, p) is the singlet and S (e, p) the triplet component, we can try solutions of a pure singlet form, Q(e, p) ^ 
while S(e, p) = 0, or of a pure triplet form, Q(e, p) = while S(s, p) 7^ 0. At the same time, the symmetry of M 
implies that 

Q(e,p) = -Q(-£,-p) and S"* (e, p) = S 4 (-£, -p) for i = 1,2,3 . (54) 
In the first case of a pure spin singlet, the self-consistency reduces to 

£> eff (e-e',p-p';[Q])~ d P' 



Q (e,P) = 3«A T > / Ei — Si— —2, (55) 

ie'-Vp +iQ(s',p') (2tt) 2 

while in the triplet case, we obtain 

, ^ f D eS (e-s',p-p';[S}) dp / 
5(e,p) = iA 2 rV / Ekt J -^— ^ (56) 

Both the equations and (|55)) have non-trivial solutions. Let us for the moment restrict ourselves to the singlet 
component. 
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First, we write the solution for Q in the form of a diagonal matrix in L-space, cf. Eq. ([16 

Q 1 (e,p) 



( £ 'P)-( ' Q 2 (e, p) ) ■ < 57 > 



L 



This implies that the Q L (e, p) are independent of each other, reflecting the independence of the hot-spot pairings in 
perpendicular directions in the Brillouin zone in Fig. 1. 

For each direction L, we seek for the singlet solution Q L (e,p) using the representation 

Q L (e, p) = a(e, p) signe 1 + iVpw (e, p) + iB L (e, p) S 3 (58) 

where 1 denotes the unity matrix. The matrix B L (e,p) is unity in all spaces except A- and particle- hole r-space, 
a(e, p) and w(s,p) are functions without any matrix structure. By construction, these two functions and the matrix 
B L (e, p) satisfy the parity relations 

o(e,p) = a(-s,-p), w (e, p) = to (-£, -p) , B (e,p) = -B (-e, -p) . (59) 
Further, we assume that the matrices B L (e,p) anticommute with the matrix A3, 

[b l (e, P ),A 3 } =0 . (60) 
The properties ((591 [60)) for the matrix B(e,p) allow us to write it in the form 

B L (s,p) = b(e,p)U L with U L = 1 ~f ) . (61) 

Herein, w L denotes a 2iV x 2A^ unitary matrix rotating particle-hole and flavour components and b (e, p) is a real even 
function, b(e, p) = £>(—£, — p). The matrix C/ L is unitary, U L+ U L = 1, and also satisfies U L U L = 1, implying for the 
matrix u the symmetry 

ti L = u L+ . (62) 

As a result, detu L = (detu L )* so that u L belongs effectively to the SU(2iV) group. For N = 1, we may apply the 
parametrisation [Eq. (7) of the Article] 

u L = ( ) with |A^| 2 + |A^| 2 = 1 . (63) 

Recalling the structure of the vectors Eq. (|15p . we have concluded in the Article that plays the role of the 
superconducting order parameter, whereas A^ is the order parameter for an insulating electron-hole pairing order. 
Equivalently, we may employ a parametrisation in terms of a 4-component unit vector, 

A L = (ReA^ImA^ReA^ImA^) 

= (cos 9 L cos x L , cos 9 L sin \ L , sin 8 L cos ip L , sin 9 L sin tp L ) . (64) 

The angles 8 L determine to what extent the ordering is superconducting and to what extent insulating. At 9 L = 0, 
only the insulating order exists in the L-dircction, while at 9 L = ir/2 the subsystem should manifest superconducting 
properties. The angles \ L an d f L arc the phases of the insulating and superconducting order parameters, respectively. 
In the linear spectrum approximation used so far, the directions of the vectors A L are arbitrary. 

Using Eqs. (|58l |61~T) . we reduce the self-consistency equation (|55|) to a system of three equations for a (e, p), w (e, p), 
and b (e, p), 

f(e,p)-e = 3A 2 TV / A* (e -/,P - pW.p') (65) 

V J f 2 ( e '' P') + ( V P') (! + w ( £ '> P')) + b 2 (s', p') (2n) 2 

b(e,p) = 3X 2 TT I D cS (e-ei,p-pi)b(e\p<) 

^Jf 2 (s\pi) + (vp>) 2 a + w(e>p>)) 2 +b 2 (e\n>)(2n) 2 ' 



f 2 (e>, p') + (vp')" (1+w (e>, p')) z + b 2 (e', p') {2ttY 
D cS (e -e',p- p') (vp') (1 + w (e', p')) dp' 



/ 2 (£', p') + (vp')' (1 + w (£', p'))' + b 2 (£', p') (27T)' 



(67) 
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where / (e, p) = e + a (e, p) signe and v stands for the Fermi velocity at a hot spot. In the approximation used below 
the solution weakly depends on the momentum and we omit in the following the momentum in the arguments writing 
/(e), w(e), and 6(e). 

Including the mean field in the fermion propagators, the polarisation bubble n(w,q), Eq. (|38[) rcnormalising the 
bosonic propagator D c g (e — e', p — p'), Eq. ([37]) . takes the form 

n M = -wtnJ: /(*)/(* + ") + * W ^ + "J - W Up + gj V2j j*p_ (68) 



/(e)/(e + c) + 6(e) 


6 (e + w) - 


- (P v i) 


((p + q) v 2 ) 


dp 


[f 2 (e) + (p Vl ) 2 + 62 (e)~ 




> (e + w) 


+ ((P 


+ q)v 2 ) 2 +62 (e)" 


(27T) 2 



The velocities vi and v 2 may be chosen to be those given by Eq. (fl4|) but the result does not depend on the choice of 
the hot-spot pair Also, anticipating the limit to be discussed below, we have neglected the function w{e). Since Vi 
and v 2 are linearly independent, the projections of p on vi and v 2 are a convenient choice for the two independent 
variables. As a result, we immediately see that II depends only on ui and, neglecting the irrelevant w 2 -term of the 
bare bosonic propagator, we find 

(u, q) = N ( 7 fi (w) + q 2 + a) . (69) 
Herein, the dynamic part is characterised by the function Q,{lj) 

f{e)f{e + oj) + b{e)b{e + w) \ 



V/ 2 (e) + 62 (e) V/ 2 (e + u) + V (e + w) 



(70) 



where the coupling constant is given by 



4A 2 

~f=-f—x ( 71 ) 
irv z smo 



and S is the angle between the two Fermi velocities in one hot-spot pair, 

v[ =1 = v sin {8/2) , v l y =1 = v cos {8/2) , (72) 

see also the inset of Fig. 1. All the approximations made throughout our calculations are justified assuming a small 
but nonzero angle, 

8 < 1 , (73) 

while N can be arbitrary, including the physical value of N = 1 . 

In the absence of the gap 6(e), the function O(cj) has the form of the well-known Landau damping 8 , 

n(u) = M . (74) 

The existence of the gap reduces the damping which in turn pushes the interaction D e g to become stronger. 

Equations (|65H67l l6"9"U7Tj) give us the closed system of equations for the physical quantities we are interested in. We 
can further simplify Eqs. (|65H67[) assuming that the important contributions in the integrals come from the region 
of momenta \p' ± | <C W | , where p',, and p' ± are the components of the vector p parallel and perpendicular to the 
Fermi surface. The inequality (|73|) guarantees that this is indeed the case. Moreover, this limits allows to neglect the 
function w;(e,p). 

Neglecting |_pj_| with respect to \p',,\ in the function D e g (e — e',p — p'), we integrate this function separately over 
p' u and the rest of the integrand over p' ± . As a result, the solutions / and 6 no longer depend on the momenta and we 
simplify Eqs. p)] 157 )1 to 

Pf pp {£>) + &2 ( £ /) v > 

*W - p-f> f (*"> , (76) 



where 



D[u) = nrr^ - (77) 
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At the critical point, a = and Eqs. (|75H77[) become universal: Introducing the energy scale 

r= (s^) 27rA2sin<5 (78) 

and dimensionless functions / = f/T, b = b/T, f = T/T, ft = ft/r, e = s/T, and Q = u/T, wc reduce Eqs. $M [73 
[76]) to the form 

= rV , 1 - (79) 
1(e) = rV- J - Hs>) ^, (80) 

V ^(e-e') V/ 2 (e') + 62 (F) 

ft («) = nfvfl - + ^ + + \ . (81) 

e\ VP W) + b 2 W) VP (e + Q) + V(e + Q)J 

We see that the typical values of the functions /(e) and 6(e) are of order T and thus important momenta p±_ in Eqs. 
(I65H67P of order T/v while the momenta p\\ entering the bosonic propagator D c g, Eq. ([69]) . are of order (T^f) 1 ^ 2 . 
This gives the estimate |p_i_|/|pj|| ~ (T/v 2 "/) 1 / 2 ~ (smS)/N <C 1, confirming in the limit ([73]) the validity of the 
approximations we employed when simplifying Eqs. (|65M66[) . Alternatively, we can write Eqs. (|79H8ip in the form of 
Eq. (6) of the Article, introducing the "angle" 6(e). If we wish to study the vicinity of the QCP on the metallic side, 
we can include a finite a > in Eqs. (|79H81[) replacing ft by ft + a/^T). 

Equations ([79] 180]) arc well defined at T = but lose their sense at finite T due to the formal divergence of the 
term with the frequency e' = e in the R.H.S. (ft (0) =0). As we have discussed in the beginning of this Section, the 
neglect of the quartic term in the Lagrangian L^, Eq. (|13[) . may be justified at finite temperatures by assuming that 
the function a (T) remains finite at the QCP but tends to zero in the limit T — > 0. 

One can understand this statement considering the first order correction to the coupling constant coming from the 
zero Matsubara frequency. In the first order, the zero frequency renormalisation takes the form 

'^' j - cT ' j iwT^- ,82) 

where c is a numerical coefficient. We can neglect this contribution only if a 3> Tg. At the same time, we should keep 
in mind that at finite temperatures, the antiferromagnetic transition in 2D is smeared because the thermal fluctuations 
destroy the antiferromagnetic order. The smearing of the transition means that effectively the "mass" a (T) cannot 
turn to zero at critical point at T ^ because the latter does not exist. In order words, considering fluctuations of the 
static component of 0, we have a cutoff at fco ~ VT. In order to avoid all these complications, we simply drop the term 
with e' = e from the sum over e' in Eqs. ([79l [80]) when studying these equations numerically. This is how our result 
in Fig. 3 has been obtained. Neglecting this term can lead to a somewhat lower mean field transition temperature 
into the pseudogap state To but its precise value is not very important for our present discussion. Anyway, at low 
temperatures this should be a good approximation. 

The numerical solution of Eqs. (|79H8ip is also represented at T = 0.001 in Fig. [I0l In a) one can see a finite 
function 6 (e) leading to a gap in the fermionic spectrum. 

At a finite a one can see in b) the Landau damping (linear dependence of ft (u>) on \u>\), whereas at a = the 
dependence of ft (Hi) on u> is quadratic. The latter is a consequence of the existence of the gap in the spectrum of the 
fcrmions. A more general picture is given by Fig. 3 of the Article. 

Actually, the ui 2 asymptotics of ft (ui) at small u leads to a logarithmic divergence of 6 (e) in the limit e — > 0, which 
is clear for Eqs. (|79I80[) . However, this divergency is not important for our consideration. Moreover, this logarithmic 
divergence may be cut off by taking into account the renormalisation of the coupling constant g and fermionic Green 
functions by non-zero Matsubara frequencies. This fact can be understood writing at low temperatures instead of Eq. 
([82]) the following renormalisation 

eg 2 f cPkcku 

g ^ g -2^J (n(a,) + fc* + fl ) a (83) 

As (cj) oc ui 2 at small u), one can understand from Eq. (|82p that the renormalisation of g is equivalent to the one 
of a classical three dimensional </> 4 -model. Scaling relations and exponents are very well known for this model. The 
bosonic propagator D (w, fc, a) can be written at a = in a scaling form 19 

D(u,k,0) cx k- 2+v a(k/uj), ?7« 0.04, (84) 
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FIG. 10: (a) The dimensionless quantities b and / as functions of the reduced frequency e at temperature T — 0.001 at the QCP 
(a — 0). Under the peak of the gap function b, the dynamic term in the fermionic propagator is linear, while the interaction 
with the Landau-damped bosons leads to the characteristic square root law only at larger frequencies e ~ 1. Inset: The same 
but far away from the criticality (a — 1). The gap b is considerably suppressed and the Landau damping is ineffective, resulting 
in the linear frequency dependence as in the case of free fermions. (b) The dynamic term Q of the boson propagator. At zero 
bosonic mass, the fermionic gap b suppresses the Landau damping and Q, oc uj 2 , while at a large mass a — 1 the fermion gap b 
is small and the [w[-law is restored. 



where a (x) is a function of x, with a finite a (0) and 

a (x) oc x 2 ~ v , x — > oo 

Substituting Eq. (|84|) into Eqs. (|75H77[) we can see that the logarithmic divergence of the solution b (e) is cut because 
the integral over e' converges. As the exponent rj is numerically small, the solution and all subsequent quantities are 
very close to those obtained for r\ = when the quartic term in the Lagrangian L^, Eq. (|13p . was neglected. 
As concerns the triplet order parameter S, Eq. (I53|) . we seek a solution S in the form 

S(s) = e(s)^ + , (85) 

where e 1 (e) = e l (— e), z + z = 1, and z = z + . This leads us to an equation for e 1 , i — 1,2, 3, 

e -( g )^A 2 V ( b ^ e - £, 'V P ' g) ^ ■ (86) 

e 2 + (vp') 2 + e 2 (e') (2tt) 2 

In Eq. (|86|) the propagator £) e ff is a 3 x 3 matrix. This equation differs from Eq. (|66|) by the absence of the factor 3 
in front of A 2 , which implies that the triplet solution is energetically less favorable. Therefore, we do not consider the 
triplet pairing anymore. 

In conclusion, the mean field treatment of the interaction between fermions and critical spin modes leads to a 
spontaneous breaking of the symmetry and formation of a gap. The emerging order is a mixture of a singlet d- 
wave like superconductivity and a particle-hole insulating order. Of course, we cannot fully trust mean field theories 
unless we have demonstrated that the mean field state corresponds to a the minimum and we have estimated the 
contributions coming from fluctuations. We will do this in the next Section, while in the remaining of this one we 
discuss the symmetries and properties of the new state. 



C. Symmetry of the order parameter 

The off-diagonal part Q s (in A-spacc) of the matrix Q, Eqs. ([58l IfJTj) . is the matrix order parameter containing 
the superconducting (A + ) and an insulating particle- hole pairing (A„), cf. Eq. (|63|) . This order arises as a result 
of the pairing of particles and/or holes from hot spots that are located opposite to each other on the Fermi surface. 
At the same time, the presence of S3 in Q g means that the order parameter changes sign when passing from a hot 
spot to its partner connected by the vector Q, see Fig. 1. Such an oscillation of the sign corresponds to a d- wave- 
like structure of the wave functions. Recalling that the order parameter Q originates from the spin singlet pairing, 
we thus conclude that the superconducting part of the order parameter describes cZ-wave superconductivity. A full 
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body of experimental methods for its identification exists and superconductivity is a phenomenon that can hardly be 
overlooked in transport or thermodynamic measurements. 

In contrast, the particle- hole part of this order is somehow "hidden". As it is singlet with respect to spin, no spin 
order appears. 

Now, we want to find what observable correlations arise due to the the component A_ of Q. 

Absence of charge density order. Let us begin this study by directly computing the local density modulation p (r) 
due to the electron-hole ordering. This quantity can be written in terms of the average 

22 , 
PW = EE / ((V^ P +Q L , Q ) cxp (iQ L r) + c.c) — ^ , (87) 



L=l a=l 



where ip Pt a are fermionic fields at equal times in the conventional formulation [p measured from the T-point], a stands 
for spin, and Qz, with L = 1,2 are the vectors connecting opposite hot spots ±p, cf. Fig. 1 (see also Ref. 14 ). The 
vectors Qi 2 ar e directed along the diagonals of the Brillouin zone. Using the fields x a and the definition of the vectors 
Eqs. (|15H17p , we rewrite Eq. ([87]) as 

2 2 

pW=EE/ <X« X«') eff (*Q*r) * = _i /" tr [<r 3 A 2 M>#> eff exp (ir 3 A 3 Q L r)] * , (88) 



|i|=l a=\ 



where (■■•) off stands for the averaging with the action S e g, Eq. (j4"5)) . and, by definition, Q/ = — Q_ 
Using the Green function G, Eq. (|44|) . obtain 



i(r) = iTY, J tr 



r 3 A 2 [ie-Vp + iQ (e) exp (ir 3 A 3 Q L r 



dp 



(89) 



T 2 

— ^ tr [r 3 A 2 C/ L E 3 cxp (iT 3 A 3 Q L r)] Y Q (T) 

L=l 



where 



Y (f) = fV ^^ _ S with 5 = / dpii (90) 

and U L is given by Eq. (fSTj) . The momentum 5 is the length of the arc on the Fermi surface covered by the gap. 

At first glance, the integral determining S diverges at large p\\ . This is a consequence of neglecting the curvature 
of the Fermi surface. In the next Subsection devoted to the effects of curvature we shall estimate the integral S as 

- 3A2 (») 



2Nv sin S 

A quick glance at the second line of Eq. (|89|) reveals however that the trace yields zero anyway due to the presence 
of the matrix S 3 . Therefore, the electron- hole pairing considered here does not lead to any local charge density 
modulation 

p(r) = (92) 

In fact, this reflects the <i-wave character of the pairing leading to an oscillation of the order parameter when moving 
along the Fermi surface. 

Absence of orbital currents. As concerns an average orbital current j, we start with a formula analogous to Eq. (|87[) . 

2 2 f dp 
j( r ) = - e EE / (( v ( p ) V'poV'p+Qi.a) exp (iQ L r) + c.c) — ^ , (93) 

L=l a =l J ( 27r > 



where v = de (p) /dp = [v x (p) , v (p)jJ is the velocity. Expressing the current in terms of the vector field ^ as 
Eq. (|88j) and then in terms of the order parameter U, we find 



iT 2 

j (r) = -— tr [t 3 V l A 2 [/ l 5] 3 exp ( l r 3 A 3 Q L r) + c.c.J Y (T) = . (94) 



L = l 



The trace in Eq. (|94[) vanishes due to the presence of the matrix A 3 in V, cf. Eq. (f25|) . The presence of the matrix 
E 3 under the trace would not be sufficient for the trace to vanish because V also contains this matrix. 
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Absence of dipole order. Now check whether an average dipolc moment d may arise. We write this quantity as 

2 2 



d(r) = -*EE [ ((^a^p+Q^cAexp^Qj 

L=l a=l J \\ V I 



r) + c.c 



dp 

(2^7 



or 



d M = T E/ 



tr 



T3A2 ( (is — Vp + iQ (e)) I exp (ir 3 A 3 Q L r) + c.c. 



dp 



2 

[ tr r 3 A 2V L (ie -Vp + iQ (e)) * exp (ir 3 A 3 Q L r) 



+ C.C. 



£ L = l 

T EE 

e L=l 



<ip 

dp 

(2^ 



(95) 



(96) 



tr 



r 3 A 2 V 



L -f ( £ ) + (yp) 2 + 6 2 ( e ) 
/ 2 (e) + (vp) 2 + 62 (e) 



exp (iT 3 A 3 Q L r) + c.c. 



dp 

(2^ 



= 



so that no average dipolc moments exist. Since no macroscopic scalar and vector order parameters can be constructed 
for the electron-hole pairing, we now try to compose a tensor order parameter. The natural choice is to consider the 
2x2 quadrupole tensor. 

Quadrupole order. In the model considered here, the 2x2 quadrupole moment tensor D (r) can be composed as 



2 2 f d 
D (r) = -e 2 EE / ((r p M P+ Q L> a) exp (*Q L r) + c.c) * , 

where the traceless 2x2 matrix operator I p equals 



(97) 



d z _ d 2 
dpj _ dp 1 



i)' 1 



(98) 



\ dpxdpy dpi dp 1 J 

and e is the electron charge. Proceeding in the same way as previously, we write 

2 



e L=l 



tr 



r 3 A 2 /p [ (ie- Vp + iQ (e) exp (ir 3 A 3 Q L r) + c.c 



dp 

(2^ 



(99) 



Now we can calculate the derivatives with respect to momentum p. It is clear that when neglecting the curvature 
only the derivatives along the vector V are not equal to zero. This allows us to write the quadrupole tensor D (r) in 
the form 



D (r) = J2 D L ( r ) with 



(100) 



L=l 



r 3 A 2 I 2V l L V 1 L - ( ) 5. 



where i,j £ {x, y}. The evaluation of the integral yields 

2 



(ie - Vp + iQ (e)) exp (ir 3 A 3 Q L r) + c.c. ) — — 

/ (27T 



dp 



e 5 sin 5 

7tT 



Y 2 (T) ^ tr 



L=l 



2V l L V 3 L - [V L \ Sa J ir 3 A 2 E 3 C/exp (ir 3 A 3 Q L r) + c.c. 



where 



3b 2 (e) 



(101) 



(102) 



£>0 



[P(e) + P{e)f 2 {f 2 (e) + bz(s 



,5/2 



(103) 
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Computing the trace in Eq. (|102[) and using Eqs. (|78I91|) . we find the quadrupolc moment tensor in the form 

D (r) = D (r) ( ° I) with (104) 



1 

D (r) = "^smT^ ^ ) E (Qir ~ XL) (105) 

L=l 



128e 2 iV 2 cos0 



This formula shows that the electron-hole pairing order corresponds to a quadrupolc moment oscillating in space in the 
directions Qi and Q2. This results in a chcqucrboard structure as illustrated in Fig. 2, which is a quadrupolc density 
wave (QDW). The vectors Qi,2 are directed along the diagonals of the Brillouin zone but their length is considerably 
smaller than that of Q = (tt, tt). The tensor D is diagonal for coordinate axes parallel to the boundaries of the 
Brillouin zone. Remarkably, the quadrupole density in Eq. (|105|) does not depend on the strength of the spin-fermion 
interaction. The quadrupole moment essentially depends on the angle 9, Eq. and vanishes at 9 = ir/2, which 

corresponds to the d-wave superconducting state. 



D. Effect of curvature 



The mean field equations (|79H81[) have been derived neglecting the curvature of the Fermi surface. As a result, the 
solution Q (e) does not depend on the position on the Fermi surface and one might come to the conclusion that the 
entire Fermi surface is covered by the gap. Of course, this is not the case and the Fermi surface is covered by the gap 
only near the hot spots where the spectrum has been linearised. 

In this Section, we estimate the length of the gapped part of the Fermi surface around a hot spot. 

Figure [TT1 illustrates the geometry for one of the four pairs of hot spots in the Brillouin zone (Fig. 1). Let us 
reconsider Eq. (|66|) focussing on this pair. The two hot spots are connected by the vector Q = (tt, tt). In terms of the 
original fcrmionic operators ipp,a with p close to the Fermi surface, we write the equation for the gap in the form 

b (e, P ) = 3A 2 V / ^( £ -e',p-p'-Q)^,pQ . (106) 

/ 2 (e',p') + (v P ') 2 + & 2 ( £ ',p') (2tt) 2 

How large is the region near the hot spots where it is justified to consider the gap b(e,p) as independent of the 
momenta p|| along the Fermi surface? To answer this question, let us assume that b(e',p') in the integrand on the 
right-hand side does not depend on p'. Furthermore, assuming that typical values of the momentum components p± 
perpendicular to the Fermi surface are much smaller than parallel components (to be verified a posteriori) we may 
neglect p± in the bosonic propagator D e g and perform the pj_-integral. Finally, we have to investigate the integral 

dp' 

D eS (e-e',p-p'-Q) — I (107) 

P± =p' ± =0 ITT 

and find out when it does not depend on p. 

Taking into account the curvature close to the hot spots we approximate the Fermi surface in their vicinity by 
circular arcs as shown in Fig. [TTJ The radius po of these arcs has to be chosen in such a way that the curvature of 
the arc equals the radius of the curvature of the actual Fermi surface at the hot spot. When integrating over p' x we 
neglect any non-linear dependence of the fermionic spectrum on this variable. This is justified in the weak coupling 
limit ([28]) . The centers of the circles for the two hot spots of the pair in Fig. [TT]are located at a distance 2Pq from 
each other. Let us choose coordinates according to Fig. [TT] with the center of the coordinate system in the center of 
the Brillouin zone. Then, we express the vector Q and the two points p and p' on the Fermi surface as 

Q = 2p a (P Q /p -sm(5/2) , 0), 

p = po (Po/po - sin(<S/2 + a) , - cos(<5/2 + a)) , 

p' = po( - P /po + sm(S/2 + a') , - cos(<ty2 + a ')). (108) 

The angle i5 is the angle between the Fermi velocities Vi and V2 introduced previously in Fig. 1. 
With these notations, we have close to the hot spots (|a| , \a\ ^ 1) 

(p-p' -Q) 2 ^pl \(a' + a cos S) 2 + a 2 sin 2 d] . (109) 
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FIG. 11: Fermi surface geometry in the vicinity of a pair of hot spots. 



Since ~ poa and p'^ ~ Poa', the integral (| 107[) takes the form 

1 dp' {l 
jtt (e - e') + a + (p 1 ,, + p\\ cosS) 2 + pf, sin 2 5 27r 
1 dp' l{ 
(e - e') + a + p',f + pf, sin 2 S 2tt 



(110) 



From this intermediate result, we understand that the curvature could be fully neglected only if 6 = [or 5 = 7r]. 
This limit corresponds to an effective nesting around the hot spots and the integration over p'^ in Eq. f| 1 10|) leads to 
Eq. (fTTj) as derived previously. For finite sin 6 corresponding to non-collinear Fermi velocities at adjacent hot spots, 
we cannot neglect the last term in the denominator of Eq. (|110p so that the solution b of Eq. (|106p does indeed depend 
on vanishing far away from the hot spots. 

The exact solution can be found only numerically but an estimate for the characteristic length of the arc under the 
gap can easily be done. This estimate follows from comparing the last term pf sin 2 5 in the denominator with the first 
two. It is clear that the solution b is independent of p» so long as this last term in the denominator can be neglected. 
Hence, we estimate the length S of the gapped region as 

5 sin 5 = 2 [( 7 ft(e - e') + a)} 1/2 . (Ill) 

As we have seen in the previous subsections, typical values of f2(e — e') are of order T, Eq. ([751) . Thus, at the critical 
point, a = 0, we come to the estimate already given in Eq. (|9ip . 



s M2 



2Nv sin 5 



We have assumed throughout the mean field scheme that the essential \p± \ are much smaller than |p|| |. As \p±\ ~ T/v 
and |p|| | ~ S we obtain 

^ Tjv = 3-7T sin 2 8 
|P||| 3 327V ' v ' 

which justifies the approximation. 

Making the parameter a finite increases the area S under the gap until the entire Fermi surface gets covered by the 
gap at sufficiently large a. For the length of the gapped region S (a) at finite a, we find 

S(a)=S(0) + ^-, (113) 
smo 

where S (0) is given by Eq. ([91j) . This is clearly not more than a qualitative estimate but it shows the increase of the 
gapped region. The value of the gap b itself decreases when increasing a, though. 

All formulas obtained in the mean field approximations could also be justified in the limit of large N for any sin 6. 
However, this limit is not useful in establishing of the validity of a u-model considered in the next Section and cannot 
substitute the inequality (f?3"l) . 

The bottom line of this Subsection is that, whenever integrals J dpu appear, they should be simply replaced by the 
value S, Eq. OT. 
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III. FREE ENERGY IN THE MEAN FIELD APPROXIMATION AND DERIVATION OF THE 

NON-LINEAR cr-MODEL 

A. Calculation of the free energy. 

All the results of the previous Section have been obtained in the mean field approximation. Now, we would like 
to calculate the energy of the ground state and consider fluctuations. This is certainly important for checking the 
applicability of the mean field scheme but, at the same time, fluctuations enrich the physics and lead to new phase 
transitions. Therefore, after having established the mean field picture, the investigation of their role is the next step 
in the study of the spin-fermion model. However, let us first concentrate on calculating the free energy 

An important peculiarity of the mean field approximation is that we not only renormalise sclf-consistently the 
fermionic Green functions, Fig. [21 but also the bosonic propagator D is renormalised by including polarisation 
bubbles, Fig. [5J which in turn contains the renormalised fermionic Green functions. In principle, we would encounter 
the pseudogap state also just using the Landau damped form 8,9 ' 10 ' 14 of the propagator D, i.e. inserting the bare 
Green functions into Eq. ([38]) . However, this approximation would be incorrect now and, even though a finite solution 
for b (e) exists in this case, the values of the gap b (e) would be considerably smaller. As we have seen, including 
the self-consistent rcnormalisation of the fermion and boson modes into the mean field scheme is not difficult but the 
calculation of the ground state energy and the derivation of a field theory for the fluctuations requires a certain care. 

The most useful feature of the diagrammatic technique is the possibility to replace parts of diagrams by "blocks" 
(see, e.g. Ref. 19 ). In other words, one can calculate vertices or self-energies within a certain approximation and 
eventually work with these quantities rather than with sums of diagram classes. While such a scheme is possible for 
the calculation of Green functions, it does not apply as directly for computation of the free energy. This problem can 
be overcome by considering a derivative of the free energy or of In Z with Z being the partition function, Eq. pip . 
Here, let us therefore write its dependence on the coupling constant A explicitly We have 

2 d\nZ x _ /S int [\|/]exp{-S [*]}£>* 

dX 2 /exp{-S [*]}£>* ■ 1 ' 

Then, the average of (Si nt entering Eq. (|114|) is calculated using the block structure of our theory. The partition 
function Z x is recovered at the end as 

A 2 A' 

^Z X = { ^-dX< 2 +\nZ° . (115) 
Jo dx 

Within the mean field approach of the previous Section, we have to insert the Green function G A , Eq. ([5"Tj) . and the 
propagator D x s , Eq. ([37l [38]) . both depending on A, into \nZ x . In the lowest order in the perturbation theory, we 
have for the derivative entering Eq. (|115|) 

d\nZ x 1 .c-^ f , r „ , . ^. „ , , _, t , dpdq 

-^ = —TY2 j £(^q)tr[Go(£,p)a%G (£ + W ,p + q)a%] ^ (116) 

where the bare A- independent bosonic and fermionic propagators D and Go have been defined in Eq. (1351 136[) . The 
self-consistency scheme is equivalent to the summation of non-intersecting diagrams for the self-energy of fermions, 
Fig. [9j and polarisation loops consisting of the dressed fermion lines for the self-energy of the bosons, Fig. [8j As a 
result of such a summation, one comes instead of Eq. (| L 16[) to 

d\nZ x VdF x l^v^ f n i , , r „, , ^ . _ t 1 dpdq 



dX 2 T d\ 2 



= -\ T Y2j D eff q) tr [G (a, p) a%G (e + w, p + q) <?%] (117) 



where the Green functions D x s (ui, q) and G (e, p) are given by Eqs. ([37l [38]) and ([5T1 [58]) . F x is the free energy per 
unit volume, and V is the volume of the system. 

Using the fact that the spin blocks of the Green functions G are proportional to unit matrices we obtain 

T dA^ = 4 J D ^( x - x ')^\^iG{X - X'^GiX 1 ^ X)]dXdX' . (118) 

Using Eqs. ([5T1 [551 l6"Tj) and subtracting the derivative of the free energy for the solution / (e) of Eq. ([79"]) with 
6(e) = 0, we obtain for the difference of the free energies AF X 

^|rP = E ( £ ( £ ) - /o 00 + e (1 - sin 9 (e))) , (119) 

£>0 
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FIG. 12: Plots for the free energy — AF, Eq. (fT20)) : (a) The dimensionless quantity —A/ = -(8NT 2 S/v)' 1 AF as a function 
of the reduced temperature T and bosonic mass a = a/{^Y), (b) —Af as a function of T at zero mass, and (c) —Af/T as a 
function of T at a — 0, showing the asymptotic behaviour —Af/T oc InT for small T. 



where £ (e) = \J f 2 (e) + b 2 (e) and sin® = f (e) /£ (e). Now, we are in the position to write down a formula for 
the difference AF between the free energy in the case of the solution with a gap b (e) ^ of Eqs. ([79l l80|) and the 
free energy for the gapless solution [6 (e) = 0]. Measuring all quantities with the dimension of energy that enter this 
formula in units of T, Eq. ([78]) . we find 

AF = -8AT 2 ^T^ [(S (e) - /„ (e) +e (l - sin 6 (e)))] (120) 

e>0 

Figure [12] displays AF as a function of the reduced temperature T and the bosonic mass a. This function is clearly 
negative, which shows that the solution with b ^ is energetically more favourable than the one without the gap 
\b(e) = 0]. Although the function AF vanishes in the limit r^O the function F/T logarithmically diverges. As the 
latter enters the weight of the states, we conclude that the pscudogap does exist below the critical temperature Tq. 



B. Non-linear a-model. 



The mean field solution Q(e), Eq. (|58l l65| l66|) . does not depend on the coordinates. Once we know Q(e), we 
can write the free energy F in this approximation. At the same time, the solution Q (e) is infinitely degenerate 
because of the arbitrariness of the SU (2A)-matrix U in Eq. (|6ip . As a result of the spontaneous breaking of the 
symmetry, gapless Goldstone modes exist in the system and strongly affect properties of the 2D model under study. 
The mean field theory developed so far does not take into account these fluctuations and is thus insufficient for 
describing important effects. Therefore, we have to go beyond the mean field theory and derive a field theory that 
finally describes the gapless excitations. 

A similar problem is encountered in the Anderson localisation theory and can be successfully treated with the help 
of a supermatrix non-linear ct- model. The mean field theory for the disorder problem corresponds to the self-consistent 
Born approximation but the latter in not sufficient for studying the localisation and this is why one should study 
fluctuations with the help of the c-model. The supersymmetry is needed in order to be able to average over the 
quenched disorder in the very beginning of the calculations before doing any approximations. 

For the spin-fermion model considered here the integration over the field <fi can be done exactly (neglecting the (f) 4 
term in L^, Eq. (fTB")) '), leading to the 'J' 4 interaction, Eq. (j3"3"|) . Following the scheme developed in the localisation 
theory, we would single out slow pairs of the the fields \Jj and decouple them by integration over an auxiliary matrix 
field Q (Hubbard-Stratonovich transformation). Then, integrating over the field <F, we would obtain a functional 
of the field Q, for which the next step would be to find the saddle point equations. These equations coincide with 
the mean field equations. In the final step, one should consider fluctuations near the saddle point and derive the 
non-linear a-model. 

In principle, we could follow this route if we did not have to take into account the gap in the spectrum of the 
fcrmionic Green functions when calculating the polarisation loops that renormalise the bosonic propagator -D e ff, Eq. 
(1371 [581 . The problem is that the gap should be found self-consistcntly, so that we cannot choose the self-consistent 
propagator D e s as the interaction in the original model. Since there is no simple way to include the self-consistent 
renormalisation of D c s into the scheme based on the Hubbard-Stratonovich transformation, we have to proceed in a 
somewhat different way. 

We assume here that the entire volume in {A}-space is subdivided into subspaces {A/} and write the free energy 
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of {X} space as 



*W = £*{*,>> ( 121 ) 



where F{ Xf } is the free energy calculated in the mean field approximation but the integration over X in this formula 
is extended over the region Xf. At the energy minimum, the gap b(e) is the same in each region but the matrices 
Iff, Eq. (|61[) . can be completely different provided the regions {Xf} are isolated from each other. An interaction 
between the regions will couple U / with one another and the free energy counted from the mean field energies will be 
proportional to the sum of differences between the matrices Uf. 

Following this idea we assume that the free energy functional F for the variation of U is given by the following 
expression 

F [u] = F [u + Su] -F[u] , (122) 

where 



Su (X) dX = . (123) 

However, in order to proceed with the derivation of F [u] we should find first the proper form of the mean field energy 
F. The derivative of this energy has been written in Eq. (| 1 18|) and we have to reconstruct the energy itself. However, 
writing the free energy in a form an integral over A 2 of the expression is not helpful because we would have to derive 
the cr-model for all A, assuming at the same time that the energy of the excitations are smaller that the gaps for all A 
including the smallest values of A. Unfortunately, writing the explicit form of the function whose derivative is given 
by Eq. (|118p is not simple. 

In order to avoid complicated derivations, we simply replace the function F x in Eq. (| 1 18(1 by another function F A , 
such that 



V d¥ x 3 



T d\ 2 4 
where the propagator B A is introduced as 



(X - X') tr [EiG (X - X') SiG (X' - X)} dXdX' , (124) 



d(X 2 D x ff ) 

dX 2 ^ 



This replacement is not justified by a small parameter but we do not expect that using Eq. (I124p instead of Eq. 
(|118p may be dangerous. First, the propagators D x s and ED A have the same asymptotics in the limit of large u,vk 
much exceeding the energy scale T, because in this limit the polarisation loops can be neglected. Moreover, they have 
the same asymptotics in the opposite limit of ui, vk much smaller than T because the function f2 (w), Eq. (|70[) . can 
be estimated in this limit as (<j) cx lj 2 /T and the propagator D x s , Eq. (|69|) . does not depend on A being equal 
according to Eq. (| 125(1 to B A . Therefore, the function D A can serve as interpolation formula for the propagator D x s . 
As the gapless fluctuations arise due to the degeneracy of the order parameter, the replacement of D x s by H> A can 
numerically change coefficients in front of the gradient terms in the cr-model but not the form of the cr-model iteself. 
The new approximate coefficients will have the same order of magnitudes as the exact ones. 
First, we guess the form of the function F writing it in the form 



— F A = 
T 12 



l - J tr [P x (X, X') \ 2 D x cS (X - X') P x (X\ X)] dXdX' - ^Tr In («„ - i\ 2 {D X S P X ) Ei) , (126) 



and the matrix P (X, X') is related to the solution Q (X, X 1 ) of the self-consistency equations ([551 [rJ51 |6"T|) as 

X 2 D x tf (X - X') P x {X, X') Ei = Q (X, X 1 ) (127) 
In Eq. (TT26| the symbol (D x g P x ) means 

(D x s P x ) - D x s (X - X') P x (X, X') (128) 



and trace Tr includes both the trace of matrices and integration over X. 
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So, the free energy F A is composed as a functional of the solution Q of the self-consistency equations and the 
propagator D A ff (X — X') containing in the polarisation loops the Green functions G, Eq. (|5"T|) . with the same Q. 
In order to check that F A is really the correct free energy we calculate its derivative 



V d¥ 
TdX 1 



G 



tr 



(P x (X, X') - 3i£iG A (X, X')) D x s (X-X 



,. dP x (X',X) 



d\ 2 



dXdX' 



(129) 



~ J tr [(P x (X, X') - 6iSiG A (X, X')) B x (X - X 1 ) P (X', X)] dXdX' 



Using Eq. (|46l l53|) one can understand easily that the integrand in the first integral in Eq. (|129|) (written for matrices 
having unity spin blocks) vanishes. Using the same self-consistency equation for calculation of the second integral in 
Eq. ([T2T)]) we come to Eq. (Tl2"3|) . 

With the help of Eqs. (TT27I \V28\i . we rewrite Eq. (fT26|) as 



T 12 J X 2 D x ff (X-X') 2 v ^' 



(130) 



Next, we calculate the energy of a slowly varying in space and time matrix U for each configuration of these matrices. 
We represent the matrix Q in the form 



Q (e, X) =a (e) signe + ib (e) U (X) , 



(131) 



with the matrix U defined in Eq. (|6T[) . 

Alternatively, the free energy functional ¥[u), Eq. (|130[) . could be obtained by decoupling the interaction term, 
Eq. (j3"3"| . singling out slow pairs of the fields and decoupling them with the help of the Hubbard-Stratonovich 
transformation. It should be kept in mind though that following such a procedure one would have to use as the 
bare interaction the propagator A 2 D A ff (X — X') with the preliminarily found gap function 6(e). Then, looking for 
the saddle point and expanding in fluctuations one should not differentiate D x s (X — X 1 ) with respect to the gap 
function. 

The first term in Eq. (|130p does not give any contribution for Q (e,X) from Eq. (|131|) because U (X) drops out. 
As concerns the second one, we write 



SQ (s,X) = ib(e)5U{X) 



(132) 



and expand the second term in SU (X) up to second order terms in space and time derivatives. We express F [u] in 
the form 



FN 



1 



tr 



SQ (e;-u,-q)G(e 



lj q 
2' P ~2 



dpdq 



(133) 



where 



G(e,p) 



*/(e)-Vp-6(e)A a 



and / (e) = (|e| + a (e)) signe. 

Using the anticommutation relations 



{SU,U} = {SU,V} = 



(134) 



(135) 



and the independence of the free energy of homogeneous in space and time rotations of the matrix U we express F [u] , 
Eq. (fT22"j) . in the form 



F[ 



tr 



f (e + %) f (e - |) + (V (p+§)) (V (p-f )) + b (e - 



(136) 



be 



(V(p+f)r + & 2 (e- 



r- 



(V(p-f)) +b*(e- 



dpdq 
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Then, we expand the integrand in u> and q up to the second order and finally obtain the free energy functional F that 
describes the gapless excitations in the form of two decoupled non-linear 2 + 1 dimensional er-modcls, 

2 

FM = (137) 
S f r T du L + (X) du L (X) T (, rN 2 du L+ du L . T , 2 du L+ du L \ i tv , ooN 



T 

with it (X) belonging to the SU (27V) group. The velocities v% and are defined as 

v i=1 = {v x v y ) , v i=2 = (t; w , -v x ) (139) 

and the superscript L in Eqs. ([137H139P specifies the two different pairs of the opposite pieces of the Fermi surface. 
We see that the fluctuations of the order parameters corresponding to these pairs of the pieces are independent. 
The dimensionless constant J\ is written using the rescaled dimensionless functions, cf. Eq. (f78|) . as 



Ji = Tj2 ~ b @ - /2 ■ (140) 

This equation is valid at any temperature for T < T . As concerns the constant Jo, we write it in the limit of low 
temperatures, 

Jo = / r z de ■ ( 141 ) 

Jo ( / -2 (e - ) + 52 (£ - ) )3/ 2 

The precise numerical values are not important for the present study because we are interested mainly in universal 
properties of the model. 

The contribution of the fluctuations to the partition function takes the form 

2 f ( F L \u L ] \ 

^fluct = n j cx p ( — ) DuL ■ ( i42 ) 

As the a- model, Eq. (|138p . is effectively three dimensional at T — 0, all contributions originating from fluctuations 
are convergent at T — 0. Moreover, they are small provided the inequality (|73|) is fulfilled. We can estimate the 
typical fluctuation of the unitary matrix u writing it as 



u = 1 + Su, (Su) = (Suy = ~Su (143) 

and expanding F L in Su up to quadratic terms. Then, we can easily calculate the average square of a small fluctuation 
in the Gaussian approximation and obtain using Eqs. ([781 191| that 



1 duidk 

(144) 



Jotu 2 + Ji {ylkl + W 2fc2) (2?r) a 



The free energy functional, Eq. (| 138[) has been derived assuming that the relevant energies are lower than T, Eq. 
(|78p. which thus is the upper limit of the integration over the energies. As a result, 



=- j =- N r 

(u-lf = (Su) 2 =Y(6u ik 5u k i)~N — -sinS . (145) 

fe=i 

Eq. (|145[) shows that the average fluctuation is small at T — as sin 5 is the small parameter of our theory. At the 
same time, the large N limit does not help because the fluctuations are not sensitive to N in this limit. 

This estimate demonstrates that the mean field theory is a good approximation in the limit (f73"[) . However, the 
situation changes at finite T because the contribution of the fluctuations diverges. We can see this from Eq. (|144[) 
replacing the integral over u> by the Matsubara sum over * ne bosonic frequencies. Then, the integral over k in 

Eq. (j!44p diverges for uj = in the infrared and more sophisticated methods are required for the investigation of the 
(T-modcl. 



27 



In order to finally confirm the applicability of the cr-model Eq. (|I38p . we estimate the fluctuations of the eigenvalues 
of the matrix order parameter B (s, X). For this purpose, we represent the order parameter as 



B(e,X) 





b(e)(l + h(e,X)) 



b(e)(l + h(e,X)) 




h = h = b? 



cf. Eq. (|6ip . The ^.-excitations are gapped as seen in the corresponding free energy functional 



_S_ 

Aitv 



E 

L=l 



tr 



K g T 2 (h L (X))' 





(146) 



Herein, the constants K g , Kq 



and K\ are of order unity. 
The estimate for the second moment of the fluctuations, [h 2 ]^, is obtained similarly to Eqs. (|1441 11451) . 



w\ = — 



K g r 2 



K co 2 + K 1 (vlkl + vlkl) (2^f 



(147) 



Even at finite temperatures, we do not need to struggle with the Matsubara sum as the estimate is given in a sufficiently 
accurate way by the integral in Eq. (|147p as long as the temperature is not close to To, where the constants K gi Kq, 
and K\ are small. Then, our estimate for the fluctuations of the eigenvalues of the matrix B is 



(B - B ) 



N 

E 

k=l 



(h lk h kl ) ~ N 



vSsinS 



sin 8, 



(148) 



which is the same order of magnitude as in Eq. (|145[) . So, the fluctuations of the modulus of the order parameter are 
small at all temperatures not close to To in the limit (|T3|) . Note that, again, the large N limit does not help. 

In order to complete the derivation of the non-linear cr-model, we now take into account the curvature term violating 
the symmetry between the superconducting and QDW states. The degeneracy between these states discussed so far 
is a consequence of the linearisation of the fermion spectrum. This degeneracy is broken as soon as the curvature 
of the Fermi surface is included. We shall see that due to the curvature effects the moduli of the order parameters 
become different favouring superconductivity. This is because the curvature of the Fermi surface affects mainly the 
particle-hole channel and while the particle-particle channel is left unchanged. Formally, this difference splits the 
values of h. However, the most essential effect of curvature is the breaking of the SU (2N) symmetry of the gapless 
excitations. This is because low energy excitations are very sensitive to violations of the symmetry. In the cr-model, 
the finite curvature of the Fermi surface gives rise to an additional term. 

In order to derive this term, we write the second term in the free energy functional in Eq. (|130|) as 



^Trln^o 



with H CU rv introduced in Eq. and expand it in H CU iv The first order of the expansion does not contribute because 
the trace tr [Ut^] vanishes as U is off-diagonal in A-space. The second order of the expansion yields an additional 
term F a in the cr-model, 



F, 



T 



1 2 

iE^E 



L=l 



tr 



r 3 [V x p]' 
2m 



(if (e)+Vp + b(e)U L (X)Y 



dp 

(2^ 



dX 



(149) 



The integral over the components of the momentum parallel to the Fermi surface decouples from the integral over the 
perpendicular one. It formally diverges but we have to restrict the region of the integration by the length S closed by 



the gap at the hot spots. The rest is straightforwardly calculated and we come to the anisotropy term F a 
with 



E 



L = l 



F L 



F L 
T 



(S/2f 



Ji 



2 

E 

L=l 



tr [u L+ (X) t 3 u l (X) r 3 ] dX , 



2Q-Km 2 v 

where the constant J\ is given by Eq. (|140|) . The full free energy functional of the cr-model takes the form 



F = To + F a 



(150) 



(151) 



The anisotropy term F a contains r 3 , thus breaks the symmetry between the superconductivity and QDW and intro- 
duces a gap in the spectrum of the QDW excitations. However, this gap is small and fluctuations remain important. 
Their study is presented in the next Section. 
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IV. RENORMALISATION GROUP EQUATIONS AND SUPERCONDUCTING PHASE TRANSITION 



As we have understood, the mean field approximation works well at T = provided the inequality (|T3[) is fulfilled. 
Due to the anisotropy term F a , Eq. (|150p . the system is at T = in the superconducting state. However, fluctuations 
arc strong at finite T and can eventually destroy the superconductivity. Therefore we consider the static limit of 
the (7-model neglecting time-dependent fluctuations and fluctuations of the modulus of the order parameter. This 
approximation is enforced by the inequality (|73|) but we believe that the cr-model can be used also for an arbitrary 
shape of the Fermi surface assuming that the coefficients of the cr-model a rcnormaliscd by those fluctuations. 
It is convenient to introduce new coordinates r, 

—x,—y) , (152) 



v x v. 



and we consider here the cr-model only for one of the values L since they behave identically after rotation of the 
coordinates by 7f/2. The free energy functional F ~F/T to be studied now takes the form 

N f 

F = — / tr [Vit+Vu + k 2 t 3 u + t 3 u] dR (153) 

where 

t = ^S f K= 2_ma)/2)\ (154) 
2Ji 5 sin d mv 

In Eq. (|154p . S (a), given by Eq. (|113[) . is the length of the Fermi surface under the gap for each hot spot near 
the critical line. We see that with increasing a, which corresponds going away from the critical line, the anisotropy 
parameter k grows making the superconductivity more favorable. The minimum of the energy is achieved for u 
anticommuting with T3, {u, r 3 } = 0. In considering the fluctuations near the minimum, it may be convenient to rotate 
matrix u by tt/2, i.e. to replace u — > T\u. Then, the second term changes sign and the minimum is reached at u = 1. 
We assume that this rotation has been done. 

The standard method of studying 2D cr-model is to use the renormalisation group (RG) technique (see, e.g. Rcfs. 
16,19 ) and, proceeding in this way, all imaginable types for the symmetries of the matrices u have been studied in the 
literature. However, for completeness, we present shortly the scheme of the derivation of the RG equations for the 
parameters of the model t and k. 

Following the RG method we represent the matrix u as 

u = uqu (155) 

where uq is a "fast" varying part and u is a "slow" one. Integrating over the fast part, we obtain a renormalised 
functional. Due to the symmetries of the cr-model, the form of Eq. (f51j) is reproduced after the renormalisation but 
the constants t and n arc modified. This allows us to write RG equations for the dependence of t and k on the 
"running" momentum cutoff. 
Proceeding in this way, we have to choose a suitable parametrisation for the unitary matrix m, 

Uq = , W = W = —W . (.150) 



w 



The explicit form of w is 



W = ( -p* -I* ) 



where s = s + is a N x N hermitian matrix and p = p is complex symmetric. 

Substituting Eq. (|155p into Eq. (|153|) and expanding u in w, we write the functional T as 

T = Fu+Fi+Fk+P, with (158) 



(Vw) 2 dR, (159) 



F x = — / tr [it+Vuo*] dR, F K = — / (u£t 3 u - t 3 ) uT 3 u + dR (160) 
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where <1> = Vuu + and J- is the slow unrenormalised cr-modcl 
Integrating over fast variables one obtains 



J exp (-Jo + T\ + F K ) Duo - Z Q (cxp [- {T\ + -7 r K )]) ^ exp (~ ^(J" x + F K ) 



2 





(161) 



where 



Z = J exp (-F Q ) Du and (...)„ = Z^ 1 J (...) exp (- Jo) Du . (162) 

The integral in Eq. (|161|) gives contributions to T, thus renormalising the latter. 

The calculation of the averages in Eq. (|161| reduces to Gaussian integrals. With the help of Wick's theorem, all the 
integrals can be reduced to the averages of the type 

(w k Mw. k ) Q = -_l_(M + itrM), (163) 



(ti> t Mi)trK t M2)) = ^7fc2 tr ( MlM2 )> ( 164 ) 



where 

M = t 2 Af ' t 2 , Mi = Mi, M 2 = M 
Evaluating the averages, we obtain renormalised values t and « 

1 1 JV + 2 f 1 A 



t t 8N J k 2 (2tt) 



W 11651 

K 2 K 2 iV + 2 , /• 1 d 2 k 



t t 4N J k 2 (2tt) 



7^7^T2 - (166) 



which leads us to the RG equations 



M_N + 2 2 din (k 2 /t) _N + 2 

d£ 16nN ' dZ ~ 8ttN ' 1 j 

where ^ = — Ink. 

At small fc the RG flow for t (£) and k (^) should stop at fc ~ k, and these parameters have to be equal to the bare 
parameters to, ko in Eq. (|153|) at £ = In A c , where A" 1 is the lowest length in the theory. Then, the solutions for the 
effective temperature t («) and the anisotropy parameter k 2 read 

/ N + 2 AcV 1 , / N + 2 A c \ , , s 

* = *o(l- W oln^j and « 2 = ^1 - —t In f j K g , (168) 

where to and Ko are the bare values of the parameters t and k in Eq. (|153|) . 

We see from Eqs. (|168p that in the absence of anisotropy between the superconductivity and QDW (kq = 0) the 
logarithm diverges and the effective temperature t goes to infinity when approaching the pole. This means that the 
long range order in the pseudogap state does not exist and correlation functions decay exponentially. A finite k can 
stop the growth of t in the region when this value is not very large, t <C N. If this happens, which corresponds 
to large k, the superconductivity is stabilised (though with a power law decay of correlation functions of the order 
parameter). Decreasing no or increasing to, one approaches the region when t becomes of the order TV and we can 
expect here a transition into the disordered pseudogap phase. The critical temperature T c of this transition has a 
completely different origin with respect to the one in the BCS theory. Its ratio to the value of the gap depends on 
many parameters and definitely is not universal. 

Remarkably, the anisotropy parameter k 2 vanishes simultaneously with the divergence of t. This means that at 
the transition not only the coherence of the superconducting state is lost but the system passes into the disordered 
mixture of the superconducting and QDW states. 

Eqs. (|167p demonstrate again that the limit — > oo is not special and the divergences in the perturbation theory 
would arise anyway without being suppressed by the parameter 1/N . An interesting possibility of the supcrconductor- 
QDW transition could exist if one could change the sign of the anisotropy parameter k 2 in Eq. (| 1 53() . Certainly, it 
is not possible to the right of the antiferromagnet-normal metal transition (where a > 0). However, the pseudogap 
cannot sharply disappear after crossing the line a = to negative values of a. We thus do not exclude the possibility 
of changing the sign of k 2 in the region of negative, not very large a. If this happened the phase diagram in Fig. 6 
would contain an interesting new phase of the QDW state with a transition into the disordered pseudogap state at a 
temperature Tqdw arL d a quantum phase transition into the superconducting state at k = 0. 
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V. SMALL PARAMETER OF THE THEORY 

We have performed the calculations exercising care in keeping under control all the approximations. Let us sum- 
marize what we have assumed and what we should assume in addition in order to justify the results obtained. 
First, we have assumed that the interaction is weak, 

A 2 < vp . (169) 

Apparently, the very existence of the spin-fermion model can be justified only in this limit. Integrating over momenta, 
wc assumed that the main contribution comes from the vicinity of the Fermi surface and this is only possible provided 
the inequality (|169|) is fulfilled. This is the usual weak coupling limit of a non-ideal Fermi gas. 

The singular form of the propagator D (w, q), Eq. (|35|) . at the QCP makes the situation more complicated and the 
inequality (|169| is not sufficient to control a perturbation theory. In order to make use of the diagrammatic technique, 
it was suggested 9,10 to introduce an artificially large number N of fermionic species. This became the basis of the 
subsequent considerations of the critical behavior near the QCP in the spin-fermion model with the conclusion that 
the earlier Hertz theory 8 properly described the system. Recently, it has been realized 12,13,14 that there is a class of 
diagrams that are, against the first intuition, not suppressed by 1/N in the large N limit. 

In our work, we confirm the conclusion of Refs. 12 > 13 > 14 that a large N limit does not help in summing ladder-like 
diagrams that can interact with each other. This is analogous to the situation arising in the theory of the Anderson 
localisation where one is able to justify the use of diagrams without intersections leading to the conventional Drudc 
transport formulas but correction require the resummation of smaller diagrams leading to effective Goldstone modes 
called "diffusons" and "cooperons" 16 . We encounter such diverging contributions due to the effect of the spontaneous 
breaking of the symmetry, which generates gapless excitations. Estimating the contribution of the massive and 
massless modes we have seen that, indeed, the large N limit does not help to keep the approximations under control. 
Fortunately, we have realised that a small parameter exists in the theory, namely 8, where S is the angle between the 
Fermi velocities of two hot spots connected by the vector Q. 

Indeed, assuming 

sin<5<l (170) 

allows us to distinguish between the scales of the momenta parallel and perpendicular to the Fermi surface when 
deriving the basic equations for the order parameter and neglect the massive modes when deriving the er-model. 
Moreover, it justifies neglecting the quantum fluctuations at zero temperature and thus confirms the validity of the 
mean field approximation at T = 0. The inequality (|170[) corresponds to a tendency to the nesting of the parts of the 
Fermi surface at the hot spots with either almost antiparallel, S = 7r/2, or almost parallel, 5 = 0, Fermi velocities vi 
and V2. 

Proceeding in this way, we used the random phase approximation to calculate the effective bosonic propagator 
D c ff, Eqs. (|37l l38|) . In doing so, we did not take into account logarithmic contributions to the Fermi velocities at 
the hot spots, however. These effects have been discussed in details in previous publications 9,10,14 but neglecting 
them has to be justified. All these contributions vanish in the limit N — > oo resulting in an Eliashberg-typc theory, 
i.e. a renormalisation of propagators while neglecting corrections to vertices. In principle, we could also assume the 
N — > oo limit. However, as the latter does not help to justify the use of the cr-model, we would have to complement 
the N — ¥ oo limit with the inequality (|170|) . Of course, this would be sufficient to prove the validity of our results but 
such a model may seem too artificial. 

Actually, we are in a position to justify our results for any N including the physically relevant value of N = 1 if we 
assume 

6 < 1 (171) 

instead of (| 1 70|) . The inequality (|170[) implies almost parallel Fermi velocities at the hot spots. The case of the almost 
antiparallel velocities does not help to keep the approximations under control. The inequality (|171[) is considerably 
stronger than the smallness in Eq. (|170p but it allows to control the approximations employed during the derivation of 
the cr-model. What remains to be done is estimating the logarithmic contributions to the vertices and Fermi velocities 
as discussed in Refs. 9 > 10 ' 14 . Explicit formulas are available in these publications and, for simplicity, wc adopt them 
for our study. 

First of all, the logarithmic contributions calculated in Refs. 9 < 10 < 14 in the absence of the pseudogap diverge in the 
limit T — > 0. As the pseudogap exists and is of the order T, Eq. (|78|) . all the logarithmic contributions arc cut off at 
low energies by T. 

The renormalisation of the velocities v x and v y has been studied 9,10,14 introducing the ratio a 



a = v y /v x 



(172) 
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An RG procedure yields the equation 

da 



d£ 2kN l + a 2 ' 



(173) 



where £ is a logarithmic running variable. Eq. f| 1 73[) has a stable fixed point a = and an unstable one at a = oo. 
The general solution of Eq. (|173|) is 

where the constant C should be determined from the condition 

a(lnA) = a . (175) 
Herein, «o is the bare value of a and A is the ultraviolet cutoff of the theory so that 

C = a T^7 lnA ' 

The flow determined by Eq. (|174j) should stop at £ = In V. 

The case of almost antiparallel Fermi velocities corresponds to small bare values ao 1 and one obtains 

Equation (|176|l shows that starting from an "incomplete" nesting with antiparallel Fermi velocities, 
(ao — (7r/2 — S) <C 1), one flows towards to a stronger one. However, one cannot move far if 

(7r/2 - 5) 2^v ln F <<1 ' (177) 

which holds unless T/A becomes exponentially small. 

The case of almost parallel Fermi velocities is characterised by large 

a ~ 5 _1 > 1 

In this limit, the solution takes the form 

^^i 1 - 2^0-^4) ' (178) 

The flow moves away from the bare value but, again, cannot go far due to the low-energy cutoff T. So, for not 
exponentially small T one can neglect the renormalisation of a in both the cases of almost parallel and antiparallel 
Fermi velocities regardless of the value of N. 

The renormalisation of vertices can be expressed in terms of the parameter a. Using results of Refs. 9 ' 10 > 14 the 
vertex r^ i ^ l + for the boson-fermion interaction in the lowest order in ln(A/r) is 

IV /A 



W = 1 + 4^ arotan UJ Ill W • (179) 

Of course, the correction is always small for TV 3> 1 but it remains small for an arbitrary N provided (t»l. We see 
from Eq. (|179|) that the nesting with almost antiparallel Fermi velocities does not help to neglect the renormalisation 
of the vertex for TV = 1 but the proximity to the nesting with the almost parallel velocities does this job very well. 

The same is true for other coupling constants of the theory. For example, the bosonic field strength renormalisation 
denoted in Ref. 14 by reads in the leading order in In (A/T) 

Z * = 1 - i^V { 'a - a + (i + ^) arCtan (s) } ■ (180) 
In the limit a>l, we may expand arctan(l/a) in l/a up to (1/a) 3 term, finding 
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Again, we recover a small pre-factor of order (aJV) 1 in front of the logarithm. 
The correction to the inverse propagator D^ 1 (w,q) has also been calculated , 
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The limit a ^> 1 allows us to expand arctan(l/a) and we see that all the pre-factors in front of the logarithms are of 
the order (aN)^ 1 , which justifies the neglect of these logarithms (In addition, we should assume that the renormalised 
g is small, g <C 1, which is the case.). 

Thus, the condition (|171|) is really sufficient for proving the applicability of our theory for any number of flavours N 
including the actual physical value N = 1. We emphasize, though, that we do not expect a qualitatively different 
behaviour of the system for an arbitrary shape of the Fermi surface including scenarios in which the inequality (|171|) 
is not fulfilled. We expect that going from small angles S to angles of arbitrary size is analogous to going from the 
weak to the strong coupling limit in the theory of (conventional) superconductivity. Although it is more difficult to 
control approximations in the latter limit, the superconducting gap is even larger. 

We have seen that the major results correlate very well with properties of the high-T c cuprates. The present 
discussion leads to the conclusion that the ultimate quantitative theory of high-T c does not have to be a strong 
coupling theory (where the fixed point goes to strong coupling and thus escapes approximate treatment) and the 
great variety of physical effects can actually be obtained from the spin-fcrmion model in the limit of weak to moderate 
coupling. 



